Theoretical search for Chevrel phase based thermoelectric materials 
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We investigate the thermoelectric properties of some semiconducting Chevrel phases. Band 
structure calculations are used to compute thermopowers and to estimate of the effects of alloying 
and disorder on carrier mobility. Alloying on the Mo site with transition metals like Re, Ru or Tc 
to reach a semiconducting composition causes large changes in the electronic structure at the Fermi 
level. Such alloys are expected to have low carrier mobilities. Filling with transition metals was also 
found to be incompatible with high thermoelectric performance based on the calculated electronic 
structures. Filling with Zn, Cu, and especially with Li was found to be favorable. The calculated 
electronic structures of these filled Chevrel phases are consistent with low scattering of carriers by 
defects associated with the filling. We expect good mobility and high thermopower in materials with 
the composition close to (LijCu^MogSeg, particularly when Li-rich, and recommend this system for 
experimental investigation. 
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Widespread use of thermoelectric materials, as an al- 
ternative technology for power generation and refriger- 
ation, remains a desirable but elusive goal. The effi- 
ciency of devices based on current state-of-the-art ma- 
terials (such as PbTe, SiGe, and Bi 2 Te3/Sb 2 Te3/Bi 2 Se3 
alloys) is low. This largely restricts them to the applica- 
tions where reliability outweighs efficiency or small device 
sizes are needed. In recent years, the center of gravity of 
the search for new thermoelectric materials has moved to- 
wards crystallographically more complicated compounds. 
This is motivated in part by the fact that low values 
of thermal conductivity, needed for thermoelectric per- 
formance, are more likely in such materials, and partly 
because such structures provide more avenues for chemi- 
cal optimization. Good examples are high figure of merit 
(ZT) skutterudite and Zn-Sb thcrmoelectrics. First prin- 
ciples calculations are proving to be a useful tool in sort- 
ing out the mechanisms by which high ZT values can 
arise in these complex materials and in exploring the ef- 
fects of various chemical modifications on the thermo- 
electric properties, lira. 

A good thermoelectric material has high thermopower, 
a low lattice contribution to the themial conductivity k 
and a sizeable electric conductivityM. The thermoelec- 
tric figure of merit is ZT = <jS 2 T/k, where a is the 
electrical conductivity, S is the Seebeck coefficient (ther- 
mopower), and k = K e \ -\- Ki a tt is the thermal conduc- 
tivity (n e i and Kiatt are the electronic and lattice con- 
tributions, respectively). Current thermoelectric mate- 
rials have ZT w 1. Large values of S are typical of 
doped semiconductors, while a is large in metals. Start- 
ing from a material with low n, the above expression 
suggests that enhanced values of ZT can be obtained 
by searching for compositions which maximize the power 
factor, aS 2 , with respect to the carrier concentration n. 
An additional requirement arises from the Wiedemann- 
Franz law, which sets a rough lower bound on the Lorentz 



number L = k/ctT [L m i n — (ir 2 kB /3e)]. A value S >160 
/iV/K is thus needed for ZT >1, even if Ki at t is negli- 
gibly small. Empirical evidence and minimum thermal 
conductivity theories imply that materials having large 
numbers of atoms in the unit cell and large atomic masses 
(soft phonons), are most likely to have the requisite low 
thermal conductivities. El Here we investigate the materi- 
als in the family of molybdenum (Mo) cluster compounds 
known as "Chevrel phases" . These are based on the bi- 
naries MogXg, where X is a chalcogen (S, Se, Te). The 
crystal structure contains large voids, which may be filled 
to yield a large variety of ternary compounds with the 
general formula M^MogXg, where M can be a simple or 
transition metal atom, or a rare-earth element. This pro- 
vides opportunities for obtaining low thermal conductiv- 
ities, analogous to skutterudites, as well as a chemical 
knob for modifying the electronic properties. These ma- 
terials fulfill the requirements for low n outlined above, 
and experimental evidence-, indicates that their thermal 
conductivity is indeed IowJj 

The band structure of the unfilled MoeXg is relatively 
simple, despite the crystallographic complexity. Follow- 
ing Ref. [l6[ the band structure can be described as de- 
riving from a collection of pseudocubic MogXg clusters. 
The chalcogens form a distorted cube, Xg, and the Mo's 
occupy the face centers of the cube. The on-site energies 
of the X-p and Mo-d orbitals are close; for instance, in 
Mo 6 Se 8 E d - E p w 0.2 - 0.4 eV (in Mo 6 S 8 it is larger, 
about 0.9 eV). These orbitals are strongly hybridized: 
the hopping amplitude t p da is of the order of 2 eV. Two 
out of five d-orbitals for each Mo atom are affected by the 
pda bonding (for the xy face they are xy and 3z 2 — r 2 ). 
The total number of affected states is 8x3 + 6x2 = 36, 
of which there are 12 non bonding p states, 12 bonding, 
and 12 antibonding pda states. In unfilled MogXg, 24 
bonding and nonbonding states (48 in both spin chan- 
nels) are all occupied, and the antibonding states are 
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much higher and are never populated in the real com- 
pounds. The remaining Mo d orbitals form 18 bands 
sitting inside the gap between the non-bonding and anti- 
bonding pda bands. The dda hopping amplitude is also 
very large, w 1.5 eV, and as a result these 18 bands also 
form bonding and antibonding combinations (there are 
no nonbonding states for d-orbitals on an octahedron). 
It can be shown that the number of bonding states is 12, 
and of antibonding states 6, so there is a gap between 
36th and 37th state in the MogXg cluster. 

The total number of electrons in the pd states in 
MogXg is 8x4 + 6x6 = 68, 4 electrons short of reaching 
the above mentioned gap. Thus, one conjectures that 
doping 4 electrons in the system will make it semicon- 
ducting, and that the states of the bottom of the gap 
would be predominantly Mo d. In particular, one expects 
the bands below the gap to be sensitive to substitution 
on the metal site, but not to filling the voids. However, 
the desired doping requires substantial filling, e.g. with 
4 monovalent atoms (Li, Cu), or 2 divalent ones (Zn) per 
formula unit. 

In our calculations, the lithium (Li) filled selenide, 
Li4_5MogSe8 is found to be particularly favorable for 
thermoelectric applications from an electronic point of 
view. Other favorable choices include Cvu-sMoeSes and 
Zn2-<sMo6Se8, while transition metal fillers or alloying on 
the Mo site are found to have too low carrier mobilities to 
be candidates for thermoelectric applications.Q We con- 
centrate on the Se based compounds, as opposed to more 
common sulfides. Se has higher mass, favorable to low 
lattice heat conductance. Furthermore, selenides tend to 
be more covalent J3 which is better from the point of view 
of carrier mobility. Experimentally there is a large ho- 
mogeneity region in the Li x MoisSe8 pseudobinary phase 
diagram near this composition. U 

We use a density-functional approach within the local- 
density approximation (LDA) for the exchange and cor- 
relation to calculate the electronic structures. In accord 
with considerations presented above, Li4MogSe8 was in- 
deed found to be a semiconductor with an LDA gap of 
about 1.4 eV. For most calculations we used the gen- 
eral potential linearized augmented planejspve (LAPW) 
methocOiid with local-orbital extensionslij to relax lin- 
earization errors and treat the Ap semi-core states of Mo 
on an equal footing with the valence states. We employ 
a well- converged basis set containing ^2000 basis func- 
tions, with muffin-tin sphere radii of 2.25 a.u. for Mo and 
Se, and of 2.0 a.u. for Li. The calculations wece based 
on the reported crystal structure of Li3.2M06Se8.tj In the 
present work, we do not take into account the trigonal 
distortion that was reported for 6 < 0.4 in Ref. 14, This 
is probably associated with partial ordering of the Li ions. 

We also considered several other dopants to the 
MogSes phase, to add four electrons per formula unit 
and make it a semiconductor. These include doping with 
zinc, copper, and some transition metals at the inter- 
stitial sites, as well as substitution of ruthenium, tech- 
netium, and rhenium on the metal (Mo) site. Zn and 



Cu fillings were found to be favorable, but the Chevrel 
phases with the needed concentration have not been syn- 
thesized experimentally. Transition metal substitutions 
were found to be unfavorable due to strong carrier scat- 
tering (as discussed above, the states below the gap are 
predominantly of Mo-Mo bonding character, thus substi- 
tution at the Mo site has a drastic effect on these states) . 
A detailed account of these studies will be published else- 
where. 



> 




FIG. 1. Band structure of the Mo^Ses 

Transport properties were obtained from the calcu- 
lated LAPW electronic structures using kinetic (Bloch- 
Boltzmann) transport theory.E3 The only further approx- 
imation was that of an isotropic and energy independent 
relaxation time for electrons. Unless there is particularly 
sharp structure of the density of states near the Fermi 
level (e.g. Pd or heavy fermion materials), this approxi- 
mation is quite good. 

The structure of Chevrel compounds can be described 
as a collection of the above mentioned pseudocubic 
MogX8 clusters wiA a small symmetry lowering rhombo- 
hedral distortionHO. Several interconnecting interstices, 
which are empty sites in the binary phases, form channels 
along the three rhombohedral directions. A large vari- 
ety of ternary compounds may be synthesized by filling 
these empty sites, e.g. Li x Mo 6 Ses, and also the proto- 
type large-cation compound PbMogSs, which is obtained 
by placing Pb into the largest interstitial site at the ori- 
gin of the rhombohedral unit cell. PbMogSs and similar 
materials have been extensively studied-due to their in- 
teresting superconducting propertiesj3'E£fli3 but semicon- 
ducting compositions have attracted less interest. Unlike 
large cations, like Pb, small cations such as Li are dis- 
tributed over 12 sites arranged as two sets of concentric 
sixfold rings surrounding the large interstice at the origin 
of the rhombohedral cellOO For intermediate concen- 
trations (1 < x < 4), the relative occupancy of these two 
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rings has beeo-a-subject of debate, at least in the case 
of the sulfide EM3 However, for x = 4 it is commonly 
accepted that three Li atoms occupy the outer ring, with 
the fourth Li atom occupying the inner ring. There is 
a tendency for Li atoms to order partially, however, this 
ordering is very little studied. 

Bonding in these compounds is understood to derive 
from covalent bonds between the Mo atoms in the octahe- 
dral cluster, and between Mo and X (Ref. [l6]). The latter 
bond has mixed ionic/covalent character. As mentioned, 
the Mo-Se bonds are less ionic than Mo-S bondsJj The 
valence bands of MoeSes shown in Fig |l| are thus com- 
posed of covalent Mo-Mo and Mo-Se bonds. The conduc- 
tion band is partially occupied (four holes per unit cell) 
at the top of the valence band manifold. As mentioned, 
the top states in this band consist primarily of bonding 
Mo-Mo states along the edges of the Mo octahedra. 



0.8 







Vlo 6 Se 8 , act 
Mo 6 Se 8 cag 


jal strucl 
3 of Li 3 2 t 


ure 

^lo 6 Se 8 
















M // 



mt r x r r x 



FIG. 2. Band structure of the Mc^Ses and of the same 
compound with the Mo and Se atoms in their position in 
Li3.2Mo6Se8 (but without Li). Note the large changes due to 
relatively small changes in the Mo-Se and Mo-Mo bond lengths, 
and overall cell geometry. 

In addition to providing the four extra electrons needed 
to dope the compound to the semiconducting regime, the 
insertion of Li in the interstitial sites leads to changes in 
the structure of the Mo-Se frame, e. g., upon saturation of 
the Mo-Mo bonds, the octahedron contracts and becomes 
more regular, and the Mo-Se distances increase. The 
corresponding change in electronic structure is shown in 
Fig||, where the bands of MogSes near the valence band 
edge are compared for two crystal structures: the ex- 
perimental structure of the binary system and the struc- 
ture of the Mo-Se frame in Lis^MoeSes (but without Li). 
Fig. H indicates that the changes in the Mo-Se frame lead 
to an enhancement of the band masses. Moreover, the va- 
lence band maxima are shifted to general k-points, which 
implies that a larger number of valleys contribute to the 
transport integrals, enhancing a without degrading S. 
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FIG. 3. Band structure of the ordered rhombohedral 
LisMoeSes compared with that of a hypothetical rhombohedral 
LieMosSes. 



As mentioned, we use the experimental structure ob- 
tained for Li3.2MogSe8 by Cava et al. in our band- 
structure calculations of the LisMogSes compound, with 
the three Li atoms placed in the outer interstitial ring, 
preserving the threefold rotation symmetry. This band- 
structure is then used to compute the transport integrals 
for the semiconducting composition, under the assump- 
tion that the system follows a rigid band behavior as 
we fill the valence bands. In order to test this assump- 
tion, we also computed the band structure of Li 6 Mo6Ses. 
In this case, we placed the three additional Li atoms in 
the inner rings, once again keeping the threefold rotation 
about the trigonal axis. The results are shown in Fig. [|. 
The primary contribution to the transport properties is 
derived from the top of the valence band, which occurs 
along the R-X line in Fig. ||. The two sets of bands are 
very similar near the top of the bands, despite the very 
strong perturbation of the self-consistent potential, upon 
the introduction of three additional Li ions and three ad- 
ditional electrons. This is key because it indicates weak 
scattering due to Li disorder and vacancies, which is a 
prerequisite for obtaining reasonable mobilities and ther- 
moelectric performance. The most pronounced difference 
between these two sets of bands near the top of the bands 
is found along the T-R directions, where the eigenvalues 
are shifted downwards in the Li 6 Mo6Ses, with respect to 
the band maxima. Nevertheless, given that the pertur- 
bation in this case is substantially stronger than in the 
composition of interest [x = 4), we see Fig. [| as an indi- 
cation that rigid band behavior between the x = 3 and 
x = 4 compositions holds to a good approximation. The 
observed rigid-band behavior also shows that the crystal 
potential felt by carriers at and near the Fermi level is 
only moderately affected by the exact position of the Li 
ions. The carrier mobility is determined by the electron 
scattering time, which in turns is defined by the change 
of the crystal potential upon moving Li atoms. Corre- 
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spondingly, when the change is small, the carrier mobil- 
ity is expected to be only weakly reduced by Coulomb 
scattering due to the disorder on the Li site. The neg- 
ative effect on the mobility will be even further reduced 
by partial ordering of Li atoms, which is known to occur 
in Li4_aMo6Se8. (A phase transition to a lower symme- 
try trigonal-phase has been reported for 5 < 0.4 in this 
compound.L-H) Furthermore, although the ordering pro- 
cesses in the mixed-filler system (Cu,Li)4_,5MogSe8 have 
hardly been studied at all, one might anticipate that, say, 
CuLi 3 Mo 6 Se 8 (this compound forms), would be better 
ordered than Li 4 MogSe8. 
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FIG. 4. Calculated thermopower of Li4_aSes for 5 = 0.05, 
0.1, and 0.15. 

The calculated thermopower as a function of temper- 
ature is shown in Fig. ^ for three different doping levels. 
All three curves reach a maximum at around 200 K. At a 
doping of about 0.1 hole/cell hole/cell (corresponding to 
a concentration of 3xl0 20 holes/cm 3 ) we obtain a max- 
imum Seebeck coefficient of ~ 150 fiV/K at T ~ 250 K. 
Without quantitative data on the lattice thermal conduc- 
tivity and carrier mobility it is not possible to determine 
ZT, but it is worth mentioning that these high values at 
high carrier concentration, if combined with reasonable 
tlattj are compatible with ZT > 1, when the material is 
optimized. 

In summary, we present first principles calculations of 
filled Chevrel phases based on MogSes. The most promis- 
ing from the point of view of potential thermoelectric 
applications is Li^^MogSes, with S ~ 0.1. This mate- 
rial can possibly have values of ZT of the order of 1 or 
larger. Among the most important questions about this 
and other filled Chevrel compounds, from the point of 
view of thermoelectric applications are: (1) whether the 
effects neglected in the current study, namely the local 
distortion of the MogSes cage near Li and partial ordering 
of Li atoms will not adversely affect the carrier mobility 
(the former effect lowers the mobility, and the latter in- 
creases it), and (2) whether the alloying with Li or Li/Cu 



would be sufficiently effective in scattering heat conduct- 
ing phonons. In this regard Cu and Li filled MogSes are 
similar electronically but presumably very different vi- 
brationally, so the Li/Cu ratio may be a useful "knob" 
for controlling the thermal conductivity. We hope that 
our calculations will encourage experimentalists to look 
closer at this promising system, (Li,Cu)4_aMogSe8. 
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